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ABSTRACT 

We have investigated the formation of dust in the ejecta of Type II supernovae (SNe), 
mostly of primordial composition, to answer the question of where are the first solid 
particles formed in the universe. However, we have also considered non-zero progeni- 
tor's metallicity values up to Z = Zq. The calculations are based on standard nucle- 
ation theory and the scheme has been first tested on the well studied case of SN1987A, 
yielding results that are in agreement with the available data. We find that: i) the first 
dust grains are predominantly made of silicates, amorphous carbon (AC), magnetite, 
and corundum; ii) the largest grains are the AC ones, with sizes around 300A , whereas 
other grain types have smaller radii, around 10-20A. The grain size distribution de- 
pends somewhat on the thermodynamics of the ejecta expansion and variations in 
the results by a factor w 2 might occur within reasonable estimates of the relevant 
parameters. Also, and for the same reason, the grain size distribution, is essentially 
unaffected by metallicity changes. The predictions on the amount of dust formed are 
very robust: for Z = 0, we find that SNe with masses in the range (12-35)M Q produce 
about O.O8Af & Md £ 0.3M Q of dust/SN. The above range increases by roughly 3 
times as the metallicity is increased to solar values. We discuss the implications and 
the cosmological consequences of the results. 



1 INTRODUCTION 



Our understanding of galaxy formation is currently mak- 
ing tremendous advances and recent investigations have also 
focused on the formation of the first luminous sources (of- 
ten referred to as PopIII objects) . Several difficult questions 
arise when one deals with these peculiarly small collapsed 
objects (for a discussion see Ferrara 2000) primarily con- 
cerning the properties of their first stars and IMF (Tegmark 
et al. 1997, Abel et al. 1998, Bromm, Coppi & Larson 2000, 
Omukai & Nishi 1999, Susa & Umemura 2000, Nakamura 
& Umemura 2000, Ripamonti et al. 2000), their response to 
the energy injection of supernovae (SN) (MacLow & Ferrara 
1999, Ciardi et al. 2000), their ability to form and preserve 
enough H2 to provide the cooling for collapse (Ciardi, Fer- 
rara & Abel 2000, Haiman, Abel & Rees 2000, Machacek, 
Bryan & Abel 2000), and their contribution to the reion- 
ization (Gnedin & Ostriker 1997, Gnedin 2000, Ciardi et 
al. 2000, Ciardi et al. 2000a) and metal enrichment (Fer- 
rara, Pettini & Shchekinov 2000) of the IGM. 

In spite of this flourishing activity little attention has 
been given to the role of dust in these early epochs. At lower 
redshifts the dramatic effects of dust have been appreciated 
when estimates of the cosmic star formation rate (SFR) were 
attempted via UV/visible surveys of distant galaxies. It was 
soon realized that approaches based on the "dropout" tech- 
nique are poorly sensitive to dust-obscured galaxies. Hence, 
the SFR deduced in this way could represent a severe under- 



estimate of the actual one, if even a rather modest amount of 
dust is present in the interstellar medium of the star forming 
galaxy. Also, some galaxies could be so heavily extinguished 
that they could be completely missed from the UV/visible 
census (Cimatti et al. 1997, Ferrara et al. 1999). 

Direct indications of the existence of dust at high red- 
shift come from the reddening of background quasars; in- 
direct evidences of dust in damped Lya systems have been 
obtained from the relative gas-phase abundances of Zn and 
Cr (Pettini et al. 1997). Fall et al. (1996) have calculated 
the cosmic infrared background from dust in damped Lya 
systems, and found good agreement with FIR background 
deduced from COBE/FIRAS data, which also seem to im- 
ply the presence of dust. Recent detection of heavy elements, 
such as carbon and silicon (Lu et al. 1998, Cowie & Songaila 
1998, Ellison et al. 2000) in very low column density Lya 
clouds (log Nh 1 < 14) at redshift z ~ 3 can potentially 
indicate that dust exists also in the Lya forest: it is quite 
natural to assume that dust grains are associated with heavy 
elements. Dust in the forest clouds would be relevant to the 
understanding of their origin and association to PopIII ob- 
jects, the heavy element enrichment pattern of intergalactic 
medium, and the thermal history of Lya clouds (Ricotti & 
Gnedin 2000, Schaye et al. 2000). 

The questions that we pose here are the following. When 
was dust first formed ? Is grain formation possible starting 
from a metal-free environment ? What are the dust prop- 



© 0000 RAS 



P. Todini & A. Ferrara 



erties and amount produced ? How are these quantities af- 
fected by metallicity changes ? 

Dwek & Scalo (1980) have shown that dust injection 
in the ISM of the Galaxy from supernovae dominates other 
sources, if indeed grains can form and survive in the ejecta. 
This has become clear after the SN1987A event, in which 
dust has been unambiguously detected (Moseley et al. 1989, 
Kosaza, Hasegawa & Nomoto 1989). Indeed, the bulk of the 
refractory elements (characterized by higher melting tem- 
peratures, such as Si, Mg, Fe, Ca, Ti, Al etc.) is injected 
into the ISM by supernovae (McKee 1989). At high red- 
shift, the contribution to dust production due to evolved 
stars (M and carbon stars, Wolf-Rayet stars, red giants and 
supergiants, novae) is even more negligible or absent. The 
reason is that the typical evolutionary timescale of these 
stars ( ^ 1 Gyr) is longer than the age of the universe, 
t H = 6.6h- 1 (l+ Z y 3/2 Gyr in a EdS cosmology, if (1+z) > 5 
(adopting h — 0.65). Thus, it seems clear that if high red- 
shift dust exists, it must have been produced by Type II 
SNe, to which we then devote the rest of this study. 



2 DUST FORMATION MODEL 

2.1 Dust nucleation and accretion 

The formation of solid materials from the gas phase can oc- 
cur only from a vapor in a supersaturated state. Because 
of the existence of a well defined condensation barrier, ex- 
pressed by a corresponding "critical cluster" size, the forma- 
tion of solid particles in a gaseous medium is described as 
a two-step process: i) the formation of critical clusters; ii) 
the growth of these clusters into macroscopic dust grains. 
The classical theory of nucleation (Feder et al. 1966, Abra- 
ham 1974) gives an expression for the nucleation current, J, 
i.e. the number of clusters of critical size formed per unit 
volume and unit time in the gas: 



J = aQ 



2cr 

TTlJll 



1/2 



C\ exp 



4/i 3 



27(lnS') 2 



(1) 



where /i — 4ira^a /ksT with a the radius of molecules 
(or atoms, depending on chemical species) in the condensed 
phase; a is the specific surface energy (corresponding to sur- 
face tension in liquids), fes is the Boltzmann constant and 
T the gas temperature; mi and c\ are the mass and the con- 
centration of the monomers in the gas phase, respectively; 
Q. — (4/3)-™^ is the volume of the single molecules in the 
condensed phase, a is the sticking coefficient and S the su- 
persaturation ratio, defined below. The subsequent growth 
of the clusters occurs by accretion and is described by: 



dv 

— - — a£lvia(t), 

at 

with the condition: 



-(0) = r* = 



2crfi 
ksTlnS ' 



(2) 



(3) 



where vi is the mean velocity of monomers, r(t) is the 
cluster radius at time t, and r* is the cluster critical ra- 
dius. Equations (|lj) and (|2|) describe nucleation and growth 
of solid particles in a gas composed of a single chemical 
species (i.e. reactions of the type Fe(gas) — > Fe(solid) or 
SiO(g) — » SiO(s)). However, there are some compounds 



(like forsterite, MgzSiOi) whose nominal molecule does not 
exist in the gas phase. These compounds form directly in 
the solid phase by means of a chemical reaction with the 
reactants in the gas phase. We need to extend the the- 
ory described above to this situation. Following Kozasa & 
Hasegawa (1987, see also Hasegawa & Kozasa 1988) we con- 
sider a vapor in a supersaturated state. In this vapor grains 
condense homologously via the reaction: 



ViAi — solid compound, 



(4) 



where Ai's represent the chemical species of reactants and 
products in the gas phase and Vi's are stoichiometric co- 
efficients, which are positive for reactants and negative for 
products respectively. We make the following assumptions: 
i) the rates of nucleation and grain growth are controlled 
by a single chemical species, referred to as a key species, ii) 
the key species corresponds to the reactant with the least 
collisional frequency onto a target cluster. In this case, eqs. 
(hi) and (0) become: 



J = a£l 



and 
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dv 

— = afivifcCifc(t), 

where mik,Cik and vik are the mass, concentration and 
mean velocity of monomers of key species, respectively. In 
this case the supersaturation ratio is expressed by: 

AG r 



In S 



,rr +E^ ln ^ 



(7) 



where Pi is the partial pressure of the i-th specie, R is the 
gas constant and AG r is the Gibbs free energy for the reac- 
tion (|). 

We investigate the formation of the following solid com- 
pounds: AI2O3 (corundum), iron, Fe^Oi (magnetite), 
MgSiOz (enstatite), Mg^SiO^ (forsterite) and amorphous 
carbon grains (ACG). These compounds are constituted by 
the most abundant heavy elements in the ejecta. Numer- 
ical constants used in our calculations are summarized in 
Tab. 1. The value of the sticking coefficient a is set equal 
to 1 for all reactions; we have checked that the final results 
are insensitive to a different choice in the plausible range 
a = 0.01 -j- 1). 



2.2 Supernova model 

We now describe the adopted model for the SN ejecta. Before 
the explosion the progenitor develops the standard "onion 
skin" stratified structure, with a hydrogen-rich envelope, a 
helium layer, and several thinner heavy element layers up to 
a Fe — Ni core. During the explosion a shock wave propa- 
gates through the layers, reheats the gas and triggers the ex- 
plosive nucleosynthesis phase. This phase lasts for few hours, 
then expansion cools the gas and the thermonuclear reac- 
tions turn off. After the explosion the SN starts to expand 
homologously, with velocity v oc R, where R is the distance 
from the center. During the first weeks Ray leigh- Taylor in- 
stabilities cause the mixing of the internal layers (Fryxell, 
Miiller & Arnett 1991). The early emergence of X-rays and 
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Table 1. Chemical reactions and numerical constants used in dust formation calculations. 
"Kozasa, Hasegawa & Nomoto 1989; Kozasa et al. 1996; c Hasegawa & Kozasa 1988. 



solid compound 


chemical reaction 


a [erg] [cm] 2 


a [10~ 8 cm] 


ACG 


C(g) -» C(s) 


1400" 


1.28 


AI2O3 


2AI + 30 -► AI2O3 


690 a 


1.72 


Fe 


Fe(g) -v Fe(s) 


1800 c 


1.41 


Fe-iOi 


3Fe + 40 -> Fe 3 4 


410 a 


1.80 


MgSi0 3 


Mg + SiO + 20 -> MgSi0 3 


400 6 


2.32 


Mg 2 Si04, 


2Mg + SiO + 30 -> Mg 2 Si0 4 


436 a 


2.05 



7-rays observed in SN 1987A (Itoh et al. 1987; Kumagai 
et al. 1988) can be explained if radioactive 56 Co is mixed 
from the internal regions of the star into the external ones; 
more precisely, observations suggest mixing of the materi- 
als in the ejecta at least up to the outer edge of the helium 
layer. Dust grains are formed by heavy elements so we fo- 
cus on the volume that containing them, i.e. the sphere of 
radius R, defined as the radius of the outer edge of the He- 
rich layer. It is thought that mixing forms clumps of heavy 
elements embedded in the He-rich layer. As a first approx- 
imation, we assume that mixing is complete, and that the 
gas has uniform density and temperature in the considered 
volume at any given time. 

Photometric observations have shown that a SN emits 
typically 10 49 erg in electromagnetic energy, but current the- 
oretical models predict kinetic energies E^ in ~ 1 jl erg. The 
expansion velocity v is then given by: v — \J ' Ekin/Mtot, 
where M tot is the total mass ejected by the SN. We take 
the chemical composition of the expelled gas from the re- 
sults of Woosley & Weaver 1995 (hereafter WW95), apart 
from the specific case of SN 1987A, see below. They de- 
termine the nucleosynthetic yields of isotopes lighter than 
A=66 (Zinc) for a grid of stellar masses and metallicities 
including stars in the mass range 11 — 40 Mq and metallici- 
ties (Z/Z e ) = 0, 1CT 4 , 0.01, 0.1, 1. They also give the values 
for Ekin and M to t for all the SN models considered. The 
range 11 — 40 Mq is the most relevant mass range for the 
production of heavy elements. In fact, stars with mass be- 
tween 8 and 11 Mq are characterized by very thin heavy 
element layers, whereas stars heavier than 40 Mq might be 
rare and give rise to a black hole partially swallowing the 
nucleosynthetic products (Maeder 1992). 

Expansion of the ejecta leads to cooling of the gas. We 
have already mentioned that the radiation losses are only a 
few percent of the total internal energy; their contribution 
to cooling is even smaller in the first week after the explosion 
due to the high opacity which prevents photons from escap- 
ing from the inner regions. Therefore, it is a good approxi- 
mation to assume that the expansion is adiabatic, although 
this hypothesis becomes less correct in the advanced evo- 
lutionary stages. Radiation losses are also partly balanced 
by the heating provided by radioactive decay (especially of 
Ni -v' 6 Co -V' 6 Fe). We neglect here these complications 
and assume that the expansion is adiabatic. In this case (for 
a perfect gas) the temperature evolution is given by 

3(1-7) 



7 is the adiabatic index, T, and Ri are the temperature and 
the radius at the beginning of the computation, t is the time 
elapsed from this initial epoch. At the temperatures of in- 
terest (T < 6000 K) the gas density is ^ 10 8 atoms cm~ 3 
so the use of the perfect gas law is well justified. We set the 
values of Ri, Ti and 7 as follows. From photometric obser- 
vations of SN 1987A (Catchpole et al. 1987) it is deduced 
that photosphere and the outer edge of Jfe-rich layer overlap 
~ 70 days after the explosion, when the photospheric tem- 
perature is 5400 K and the radius R = 1.6 x 10 lj cm. For the 
adiabatic index we take 7 = 1.25 as in Kozasa, Hasegawa 
& Nomoto (1989). We generally use these fiducial values in 
our models, but we will consider the effects of varying the 
values of Ri and 7 when discussing the results. 



2.3 Molecule formation 

In the ejecta of SN 1987A molecules of CO were detected for 
the first time in a SN (Meikle et al. 1989; Meikle et al. 1993; 
Bouchet & Danziger 1993) together with SiO (Aitken et 
al. 1988; Bouchet et al. 1991). These two molecules are very 
important for our study because carbon atoms bound in 
CO are not available to form ACG, whereas SiO molecules 
take part in the chemical path leading to the formation 
of MgSiOs and Mg2SiOi. We investigate the process of 
molecule formation under the assumption of chemical equi- 
librium. In the absence of grains, molecular formation in the 
gas phase must have been initiated by radiative processes. 
We assume that formation of CO is dominated by radiative 
association (Lepp, Dalgarno & McCray 1990; Liu, Dalgarno 
&Lepp 1992): 



C + O 



CO + hv, 



with a rate coefficient K ra (CO). The main destruction pro- 
cess of CO is the impact with the energetic electrons pro- 
duced by the radioactive decay of 56 Co (Liu & Dalgarno 
1995): 



CO + e 



C + O + e, 



with a rate coefficient K r d(CO). In steady state, the abun- 
dance of CO is given by: 



n (CO) = ^°S »(C)n(Q). 



(8) 



T = T t [ 1 + — t 
Hi 



K rd {CO) 

Analogously, radiative association is the most important 
mechanism for the formation of SiO (Liu & Dalgarno 1996) 
via the reaction: 



Si + O 



SiO + hv, 
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with rate coefficient K ra (SiO). Silicon monoxide is mainly 
destroyed by impact with energetic electrons: 



SiO + e 



Si + + e, 



with rate coefficient K r d(SiO) and by charge transfer with 
positive ions of Ne (Liu & Dalgarno 1996): 



SiO + Ne 1 



Si + + + Ne, 



with rate coefficient K c t{SiO). The abundance of SiO is 
given by: 



n(SiO) 



K ra (SiO) 



-n(Si)n(0). (9) 



K rd (SiO) + K ct {SiO)n{Ne+) 
We assume for the rate coefficients the following values: 
K ra {CO) = 1(T 16 x (-0.0398 + 1.25T 4 - 1.46T 4 2 + 



0.88T 4 3 



O.2IT4) cmV 



(Gearhart, Wheeler & Swartz 1999, Dalgarno, Du & You 
1990), where T 4 = T/10 4 K. 

K ra {SiO) = 5.52 x i - 18 T 031 cmV 1 , 

(Andreazza et al. 1995, Liu & Dalgarno 1996). 

K ct (SiO) = 2 x 10~ 12 cmV 1 , 

(Liu & Dalgarno 1996). The ejecta are only moderately ion- 



10" 



hence, we assume 



ized with fractional ionization 
n{Ne+) =0.02n(7Ve). 

To calculate the rate coefficient K r d(CO), we assume the 
same rate of destruction by energetic electron impact for 
CO as for SiO, i.e. K rd {CO) = K rd (SiO). High energy X- 
rays and 7-rays produced by the chain of radioactive decay 
H^i — > 27 C° — * 26 Fe interact by Compton scattering with 
the electrons in the ejecta. The average energy deposition 
rate per particle in the ejecta is (Woosley, Pinto & Hart- 
mann 1989): 

56 Co) 



L 7 = 7.5 x 10 



.JVi( 



Ntot 



-( J E 7 )/ 7 (^ 56 ) 



exp 



(. MeV s" 

T56 



Ni( 56 Co) is the total number of atoms of 56 Co in the 
ejecta; N to t is the total number of particles in the gas; 
(Ej) = 3.57 MeV is the mean energy of 7-ray released by 
each decay; T56 = 111.26 days is the e-folding time of 56 Co. 
The deposition function, / 7 is proportional to the fraction 
of trapped 7-photons, 

MKso) = 1 - eM-KaaMto/t) 2 }, 

with <j> = 0(t o ), the mass column density of the ejecta at 
some fiducial time t , and Kse an average opacity for 5 Co 
decay 7-rays. We assume 4> — 7 x 10 4 g cm -2 at t = 10 6 
s and 7^56 = 0.033 cm 2 g _1 . These values are appropriate 
for SN 1987A but, lacking more detailed information, we 
extrapolate them to all our models. Finally, the estimated 
destruction rate of CO and SiO by energetic electron impact 
is 

K rd {CO) = K rd {SiO) = h. s -\ 

Wd 
with Wd being the mean energy per dissociation, defined as 
the energy of primary electrons divided by the number of 



molecule dissociations (Liu & Victor 1994). For a fractional 
ionization of the gas « 10~ 2 , Wd = 152 eV. 



3 A TEST CASE: SN 1987A 

To test and calibrate our dust formation model, we first ap- 
ply it to SN 1987A, a case for which firm evidences of dust 
formation have been collected. In view of this test, we briefly 
summarize the relevant observational results. 



3.1 Observational results 

The most relevant evidence of newly formed dust grains in 
the ejecta of SN 1987A is the blue shift of the line pro- 
files. Spectroscopic observations detected this change be- 
tween August 1988 and March 1989 (Lucy et al. 1989, Lucy 
et al. 1991). This effect is likely to be caused by the larger at- 
tenuation suffered by radiation received from receding mat- 
ter due to dust grains distributed in the ejecta. The conden- 
sation efficiency, i.e. the dust mass expressed as a fraction 
of the maximum value permitted by the elemental abun- 
dances in the ejecta, derived by authors of observations from 
their data is < 10 -3 . A further interesting observation is the 
stronger fading of the [5 , J/]A1.65pim line flux relative to the 
continuum after day 530. This can be interpreted as de- 
pletion of Si from the gas phase as a consequence of the 
formation of silicate grains. If the line fading is due solely to 
depletion, than the condensation efficiency rises to > 50%. 
Another evidence of the formation of dust is an IR con- 
tinuum excess over that expected from a Planck spectrum 
fitted to the SN emission at optical wavelengths (Roche et 
al. 1989, Wooden et al. . 1993). Dust grains extinguish UV- 
visible radiation from the central energy source, re-emitting 
in the IR bands. This process might be responsible for the 
observed increase in the 10 and 20/^m fluxes around day 
350^450 (Roche et al. 1989, Bouchet & Danziger 1993, 
Meikle et al. 1993). 



3.2 Model Results 

We now turn to the main results from our nucleation numer- 
ical computations obtained by solving the above equations. 
The chemical composition of the SN 1987A ejecta is taken 
from Nomoto et al. 1991. The value of expansion velocity 
of gas is set to V — 2100 km s" 1 . This is the minimum 
expansion velocity determined from the HI Pq absorption 
trough (McGregor 1988, Nulsen et al. 1990) and is thought 
to represent the expansion velocity of the inner edge of the 
hydrogen envelope. Fig. hi shows the mass of dust formed in 
the ejecta as a function of the time elapsed since explosion 
for the different solid compounds found to be present. We 
note that there are two episodes of dust formation. In the 
first one, ACG (formation time t = 380 days) and AI2O3 
(t — 430 days) grains are formed. This process is likely to 
be responsible for the IR-excess observed at that time. In 
the second episode magnetite and silicate grains are formed, 
at about t = 600 days. Two points are worth noticing to 
this concern. The first one is that this epoch corresponds to 
the formation of the predominant fraction of dust mass of 
the SN (0.57 Mq corresponding to about 84% of the total 
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Figure 1. Dust formation as a function of time elapsed since 
explosion in SN1987A 



amount); this dust formation episode might be responsible 
for the blueshift of line profiles, observed only after day 530, 
when ACG and AI2O3 formed. 

The second point is that the formation of MgSiOz and 
MgiSiOn, might be related to the fading of the Si and Mg 
lines. The hypothesis that these elements form enstatite and 
forsterite would also be suggested by the behavior of the SiO 
molecule. The silicon monoxide rotovibrational line (v — 
1 — ► 0, A = 7.8/im) emission was detected after 160 days 
(Aitken et al. 1988) and it remained clearly detectable until 
519 days (Bouchet et al. 1991) but is no longer detected 
at 578 days (Roche et al. 1989). The time behavior of SiO 
emission can be understood by inspecting Fig. H, where we 
show the predicted SiO mass vs. the observed mass as a 
function of time. Note the rapid fall at 660 days, 70 days 
after the beginning of silicate formation, probably due to 
the depletion of Si atoms; also our model seems to slightly 
overpredict (by a factor 3.1) the amount of SiO produced at 
early times. This could be due to our simplified treatment 
of the chemical network for this molecule or inaccuracies in 
the rate coefficients. The dust formation efficiency deduced 
here is 95%, consistent with that deduced from the Si line 
fading. 

All together, we look at the above results as a satisfac- 
tory success of our model in reproducing, at least qualita- 
tively, the principal features of dust formation in SN 1987A. 



4 DUST IN PRIMORDIAL SUPERNOVAE 

In this Section we present the general results concerning dust 
formation in primordial SNe. We take the chemical compo- 
sitions of the gas from Tables 16A and 16B of WW95, and 
the relevant results are reported in Tab. 2. There the chemi- 
cal composition of the ejecta is given at 2.5 x 10 4 s after the 



Figure 2. Evolution of the SiO mass (solid line) as a function of 
time compared with the observational data (points) taken from 
Liu & Dalgarno (1996). 



explosion, when strong and electromagnetic reactions have 
ceased, but many nuclei have not yet decayed in their most 
stable form. Because dust formation occurs at 300-7-600 days 
after explosion, it is necessary to take into account the ra- 
dioactive decay of such nuclei. In the SN models of WW95 
the energy of explosion can be adjusted to give the desired 
kinetic energy of the ejecta, typically 10 erg. Following 
WW95, we explore the effects of Ekin variation by consid- 
ering a low {Ekin = 1.2 x 10 51 erg, Case A) and a high 
{Ekin — 1.9 x 10 51 erg, Case B) value for this quantity. We 
discuss separately the two cases in the following. 



4.1 Low Kinetic Energy (Case A) 

As the kinetic energy of the model is relatively low, this is 
not sufficient to completely expel the heavy elements exter- 
nal to the Ni — Fe core of the most massive SNe, and a 
variable amount of material falls back onto the core, prob- 
ably forming a neutron star or a black hole. The fallback 
will mostly affect the inner layers, containing the heaviest 
elements; as a result, progenitors with masses larger than 
« 20A/q will be prevented from forming dust. For essentially 
the same reason, above M ~ 15M©, only AGC grains are 
formed. Figs, ta-kl show the amount of dust formed as a func- 
tion of progenitor mass, and the grain composition. ACG 
are typically the first solid particles to condense, depend- 
ing on the models. The formation of these grains is quite 
fast with respect to the cooling time scale of the ejecta: 
most of the ACG dust mass forms in a narrow range of 
30 -T- 40 K around T = 1800 K. Subsequently, at a temper- 
ature of « 1600 K AI2O3 starts to condense, followed by 
Fe 3 4 , MgSiOz and Mg 2 Si0 4 at T w 1100 K. Clearly this 
sequence is governed by the condensation temperature of a 
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Table 2. Adopted chemical composition of the supernova ejecta (2.5 X 10 4 
WW95 and corrected to take into account radioactive decay. 



after explosion) for metallicity Z = 0. Data are taken from 





12M 


13M 


15M 


18M 


2OM 


22M 


25AM 
25BM 


3OAM 
30BM Q 


35AM 
35BM 


4OAM 
40BM e 


He 


4.08E+00 


4.42E+00 


4.90E+00 


5.70E+00 


6.32E+00 


7.10E+00 


7.33E+00 
7.82E+00 


9.11E+00 
9.30E+00 


8.17E+00 
1.06E+01 


8.63E+00 
1.20E+01 


C 


4.30E-02 


6.84E-02 


1.45E-01 


1.10E-01 


8.98E-02 


2.77E-01 


8.31E-02 
4.23E-01 


1.09E-01 
3.48E-01 


1.25E-09 
3.49E-01 


7.73E-10 
2.85E-01 


N 


5.41E-06 


1.00E-05 


2.46E-05 


2.13E-06 


1.76E-06 


8.39E-05 


2.39E-04 
3.29E-04 


3.64E-05 
2.74E-03 


2.13E-08 
6.29E-05 


2.66E-08 
6.17E-06 


O 


6.67E-02 


1.37E-01 


4.00E-01 


3.00E-02 


9.21E-03 


1.85E+00 


9.67E-03 
2.33E+00 


2.65E-02 
4.35E+00 


3.08E-10 
1.92E+00 


4.19E-10 
5.78E-01 


Ne 


2.61E-03 


3.04E-02 


9.28E-02 


7.72E-04 


2.12E-04 


5.57E-01 


5.79E-07 
5.10E-01 


4.66E-06 
1.04E+00 


1.37E-12 
3.76E-01 


1.95E-12 
1.10E-02 


Na 


9.83E-06 


1.90E-04 


5.68E-04 


9.57E-10 


1.01E-09 


3.90E-03 


1.04E-08 
2.85E-03 


5.04E-08 
3.88E-03 


4.02E-13 
1.86E-03 


6.37E-13 
1.07E-07 


Mg 


3.88E-03 


1.50E-02 


3.41E-02 


4.80E-06 


2.04E-06 


9.52E-02 


5.61E-08 
9.38E-02 


5.90E-07 
2.49E-01 


2.30E-13 
4.83E-02 


3.33E-13 

2.89E-04 


Al 


7.19E-05 


5.55E-04 


1.31E-03 


1.47E-10 


1.09E-10 


2.43E-03 


8.73E-10 
1.88E-03 


1.18E-08 
5.88E-03 


9.62E-14 
8.74E-04 


1.43E-13 
1.46E-08 


Si 


2.51E-02 


3.41E-02 


6.89E-02 


2.79E-09 


1.97E-09 


1.66E-01 


2.73E-09 
2.32E-01 


3.48E-08 
1.24E-01 


1.03E-09 

1.77E-03 


1.05E-10 
2.89E-08 


S 


8.66E-03 


1.35E-02 


2.61E-02 


1.18E-10 


8.87E-11 


1.09E-01 


2.07E-09 
1.06E-01 


6.36E-08 
4.59E-02 


3.58E-13 
2.82E-06 


4.22E-13 
4.32E-10 


Ca 


1.42E-03 


3.01E-03 


4.19E-03 


8.57E-13 


7.15E-13 


2.38E-02 


3.59E-10 
1.75E-02 


5.76E-09 
8.78E-03 


1.27E-26 
2.06E-09 


1.17E-17 
5.71E-11 


Ti 


1.40E-04 


1.71E-04 


1.90E-04 


3.27E-13 


6.52E-14 


3.25E-04 


1.55E-10 
3.47E-04 


3.40E-09 
4.80E-04 


2.15E-26 
1.79E-09 


1.33E-26 
1.42E-10 


Fe 


8.45E-02 


2.00E-01 


1.67E-01 


4.93E-15 


1.02E-17 


1.74E-01 


7.24E-11 
2.99E-01 


1.32E-09 
3.35E-01 


6.56E-36 
5.36E-10 


1.16E-35 
7.01E-11 


56 Co 


8.15E-02 


1.92E-01 


1.62E-01 


4.35E-20 


3.03E-23 


1.68E-01 


2.14E-16 
2.91E-01 


3.69E-15 
3.25E-01 


8.19E-46 
1.54E-15 


9.29E-46 
9.44E-16 



given material. Also, the more refractory materials preferen- 
tially end up into larger grains: this is because an earlier for- 
mation can exploit a higher concentration of the key species 
favoring the accretion process. The typical size of ACG dust 
grains is a — 300 A, whereas FaOi grains have typically 
a — 20 A and AI2O3, MgSiC>3 and Mg2SiO$, grains are 
even smaller (rs 10 A). In spite of the high condensation 
temperature, AI2O3 grains do not grow to sizes compara- 
ble to those of ACG, as their growth is limited by the low 
abundance of Al. In Figure H grain sizes of the most abun- 
dant compounds that form in supernova ejecta are shown 
for four values of Z (non zero metallicity cases are discussed 
later on). The two silicates (enstatite and forsterite) start 
to condense almost simultaneously; however, Mg2SiOi en- 
ters the supersaturation regime earlier than MgSiOz,. For 
this reason, forsterite grains grow quickly, strongly deplet- 
ing the Si (or Mg) available. Figs. 0-0 clearly show that 
even starting from a primordial composition, early SNe can 
contribute a significant amount of dust: for Case A, about 
O.O8M < M d ~ O.3M of dust/SN are produced. Interme- 
diate mass progenitors are the most efficient sources, being 
able to convert up to 2% of their mass into solid particles. 
The reason is that they synthesize a considerable amount of 
heavy elements without suffering too much from the fallback 
process mentioned above. 



4.2 High Kinetic Energy (Case B) 

In this case the kinetic energy of the explosion for the pro- 
genitor masses 25, 30, 35 and 4OM is chosen equal to 
1.9 x 10 51 erg; for lower masses the energy of explosion is the 
same as in case A. This energy is sufficient to eject also the 
inner layers, which now can provide the elements to form 
grains of various chemical composition, i.e. not only ACG 
as in case A. In Fig. pi we show the dust mass yield for Case 
B. Now SN up to masses w 35M are able to form dust. 
In addition, a SN of 3OM is able to produce about 1.3M 
of dust (4.3% of its mass). The formation sequence and the 
grain size distribution are very similar to those discussed for 
Case A. 



4.3 Effects of Ri and 7 Variations 

Up to now we have used the values of Ri and 7 deduced from 
observations on SN 1987A. These values represent a reason- 
able approximation but they might well depend on the spe- 
cific properties of the SN under examination. Therefore, as a 
sanity check, we investigate the dependence of our results on 
different choices for these paramenters. As a benchmark, we 
focus on the M = 22M SN model (Case A) and increase or 
decrease the standard value of Ri — 1.6 x 10 cm by a fac- 
tor 2.15; this corresponds to a variation of about 100 times 
in the initial volume of the ejecta. In Fig. M we compare the 
dust formation evolution and the final size of FezOi grains 
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Figure 3. Dust mass formed as a function of the SN mass (in 
the range 12M Q < M < 2OM ) for initial metallicity Z = and 
kinetic energy of the explosion E^ in = 1.2 X 10 51 erg (Case A). 
Also shown is the grain composition. 



for three values of Ri. The final masses of ACG, AI2O3 and 
Fe^Oi grains are almost unchanged in the three cases be- 
cause the gas density remains high enough for the collisional 
time scale (regulating the formation/accretion processes) of 
these materials to remain shorter than the expansion one. 
However, the behavior of silicates depends on the choice of 
Ri. As a general rule, Mg^SiO^ grains form first and grow 
faster than MgSiO*. ones, thus using up the condensable 
materials efficiently and ending up with a larger final total 
mass. However, for larger values of Ri (i.e. larger volume, 
lower gas density) this process is limited by the fact that 
the collisional scale becomes longer, thus stopping the ac- 
cretion at earlier times. The grain size distribution shifts by 
about a factor 2 as Ri is varied by a factor (2.15). Thus 
the determination of the grain size distribution is relatively 
uncertain. 

The adiabatic index 7 gives a measure of the ability of 
the gas to cool: 7 greater than the standard value of 1.25 
cause the gas to reach the supersaturation state when the 
volume of the ejecta is smaller. Therefore, a larger 7 case 
gives results similar to those obtained for the low Ri case 
discussed above. Figure H shows the dust formation evolu- 
tion and final grain size distribution of Fe^Oi grains for the 
M - 22M model (Case A) with 7 = 1.4. The similarity 
with the previous case with Ri = 7.4 x 10 14 cm and 7 = 1.25 
is evident. It has to be noted though that the variation range 
for 7 (15%) is smaller than that for Ri (50%), which might 
indicate that the dust formation process is more sensitive to 
changes in the adiabatic index than in the initial radius. 



Figure 4. Same as Fig. pi but for the SN mass range 22Mq < 
M < 3OM . 



5 EXTENSION TO HIGHER METALLICITIES 

We finally extend our results to non-primordial compositions 
by exploring the results for additional three metallicity val- 
ues Z/Zq — 10 ,10 -a ,l. The chemical composition for 
these models is also taken from WW95. We start by analyz- 
ing the dependence of the grain size distribution on metal- 
licity. This is shown in the four panels of Fig. |6|. Somewhat 
surprisingly, the dependence is almost absent, with grain 
radii ranging from 5A to 0.1 fixn for all values of Z. Also, 
the same material segregation is seen, with smaller grains 
being predominantly constituted by silicate and magnetite 
and the larger ones made by amorphous carbon. This be- 
havior can be explained by the fact that the final grain size 
is governed by the thermodynamics of the ejecta expansion 
(and therefore sensitive to Ri and 7, as already pointed out 
before), but poorly affected by the ejecta composition. 

The latter, instead, plays a more important role in de- 
termining the total amount of dust formed, as seen in Figs. 
|9l-llCJ for Case A and Case B, respectively. Moving from 
Z = to higher metallicities we observe that a large number 
of SNc contribute to dust production: a clear example of this 
is the behavior of SNe with mass 18-20 M , which increase 
their dust yield from < O.lA/ up to O.5-O.6A/ for Z — 1. 
The enhancement of dust formation is due to the fact that 
the density of heavy elements in the ejecta becomes large 
enough to allow the state of supersaturation to be reached 
more easily. This trend results in a steady increase of the 
total amount of the dust produced in the four cases (ob- 
viously, when appyling these results one has to weigh over 
the appropriate IMF); they are (2.06, 3.6, 4.5, 5.9) M for 
Z/Zq = (0, 10 _4 ,10~ 2 ,1), respectively, for Case A. How- 
ever, it is interesting to note that at each metallicity, the 
maximum amount of dust produced by a single SN varies 
little, and it is never higher than approximately lAf . Fi- 
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Figure 5. Same as in Fig. M, but for kinetic energy of the ex- 
plosion E^ in = 1.9 X 10 erg (Case B); the SN mass range is 
22Mq < M < 35Mq. For lower masses Case B gives the same 
results as Case A. 





1 
-2 
-3 


-1 
-2 
-3 


-1 
-2 
-3 

20 





I 


(a) 


1 ' ' 

t'e 3 4 ^ 






If 
lj 
If 

1 




JJ,7- 




* 


i V ^ 






(b) 




Fe 3 4 J 




1" .- 

f- 






(c) 


*'\ 


r c 3 4 J 


- 


i 


- 


i .— r"- 

1 LP 

i- i'l i 




• 


,.l , 7 



00 1500 1000 500 10 100 

Temperature [K] Grain size [A] 



10 


5 

10 



Figure 7. Dust formation (left, line types as in Fig. hi) and Fe^Oi 
grain size distribution (right) for the M = 22Mq SN model; the 
three cases refer to different values of Ri\ (a) 7.5 X 10 14 cm, (b) 
1.6 X 10 15 cm (standard value), (c) 3.4 X 10 15 cm. 
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Figure 6. Grain size distribution of ACG (open squares), Fe^Oi 
(solid triangles) and Mg2SiOi (open circles) grains for the M = 
22Mq SN model; results are give for four different metallicities 
of the progenitor. 



nally, the differences between Casa A and Case B are found 
to be minor. 



6 SUMMARY AND DISCUSSION 

We have investigated the formation of dust in Type II su- 
pernovae mostly with primordial abundances, a property 
characterizing these events in the early universe; however, 
we have also considered non-zero metallicity values up to 
Z — Zq. The calculations are based on standard nucleation 
theory and the scheme has been first tested on the well stud- 
ied case of SN1987A, yielding results that are in satisfactory 
agreement with the available data (see Section 3). The main 
results of the paper are the following: 

• The first solid particles in the universe are formed by 
Type II SNe. The dust grains are made of silicates (predom- 
inantly MgzSiOi), amorphous carbon (ACG), magnetite 
(-Fe304), and corundum (^203) and form about 300-660 
days after explosion. 

• The largest grains are the ACG, with sizes around 
300 A , whereas other grain types have smaller radii, around 
10-20 A. The grain size distribution depends considerably on 
the thermodynamics of the ejecta expansion (characterized 
by their initial radius and adiabatic index) and variations 
in the results by a factor ~ 2 might occur within the esti- 
mated range of Ri and 7. Also, and for the same reason, the 
grain size distribution, is essentially unaffected by metallic- 
ity changes. 

• The amount of dust formed is instead very robust 
to variations in Ri and 7. For Z = 0, we find that 
SN with masses in the range ( 12-35) Mq produce about 
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Figure 8. Dust formation (top) and Fe304 grain size distribu- 
tion (bottom) for the M = 22M SN with 7 = 1.4. 



Figure 10. Same of Figure JA but for Case B. 
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Figure 9. Total dust mass produced as a function of SN mass 
and different metallicity of the progenitor (Case A). 



0.08M Q < M d < 0.3M Q of dust/SN in the low kinetic en- 
ergy explosion case; slightly higher final yields are obtained 
in the high kinetic energy case. The above range increases 
by roughly 3 times as the metallicity is increased to solar 
values. 



The previous results clearly show that it is likely that 
dust has been present in the universe immediately after the 
first stars appeared. This has a large number of consequences 
that it will be necessary to study in detail. Among possi- 
ble effects, the most outstanding ones concern the opacity 
of the universe at high z, spectral distorsions in the Cos- 
mic Microwave Radiation caused by dust re-emission of ab- 
sorbed UV-optical light, catalyzation of H2 molecular hy- 
drogen formation and heavy elements depletion in the inter- 
stellar medium of pristine galaxies and in the intergalactic 
medium. 

The first two issued were already discussed in detail by 
Loeb & Haiman (1997) and Ferrara et al. (1999), and we de- 
fer the interested readers to those works for details. In short, 
the expected IGM opacity contributed by dust around the 
observed wavelength A ~ 1/im is ~ 0.13 and it rapidly in- 
creases to fa 0.35 at z — 20. The expected CMB spectral 
distortions due to high-z dust is only ~ 1.25 — 10 times 
smaller than the current COBE upper limit, but these num- 
bers might depend crucially on the formation epoch and 
abundance of dust. 

In addition to the above effects, Type II SNe can also 
initiate molecular hydrogen formation on dust grain surfaces 
rather than in the gas phase, the second process being the 
only viable in a dust free environment. As already men- 
tioned, at high redshift Type II SNe are the only possible 
sources of dust, due to the short age of the universe and 
the long evolutionary timescales characterizing more con- 
ventional dust sources, as for example evolved stars. Thanks 
to the above results we can now answer the following ques- 
tion: what is the minimum amount of dust required in order 
for the molecular hydrogen formation on grains to become 
competitive with the gas phase one ? An order-of-magnitude 
answer can be obtained by comparing the two formation 
rates. At the low densities relevant here, H2 is formed in the 
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Figure 11. Comparison between molecular hydrogen formation 
rates in the gas phase and on dust grain surfaces as a function of 
the dust-to-gas ratio and gas temperature. Above the two curves, 
indicating the equality between the rates for two different values 
of the electron fraction x e = 10~ 3 — 10~ 4 , H2 formation on grains 
dominates. 



The fate of the dust that we predict from Type II SNe 
has yet to be determined. What fraction of the grains will 
be able to survive the passage through the reverse shocks at 
which the ejecta are thermalized ? Grains in a hot gas are 
essentially destroyed via thermal sputtering, i.e. collisions 
with ions or electrons with Maxwellian velocity distribution. 
However, in spite of the still poorly understood underlying 
physics, it seems unlikely that the efficiency of grain destruc- 
tion in an adiabatic shock can be higher that 10% (Mc- 
Kee 1989). Grains are more likely to be destroyed behind 
radiative shocks, by the combined effects of a greatly en- 
hanced gas density and betatron acceleration that increases 
the grain Larmor frequency. However, if, as expected, the 
magnetic field is weak at high z, the efficiency cannot be 
very high. The grains then will follow the fate of the gas 
and will likely be expelled in the IGM, as PopIII object suf- 
fer complete blowaway of their gas (Ciardi et al. 2000). Once 
in the intergalactic space, dust might have strong influence, 
for example, on the determination of cosmological parame- 
ters via the observation of high z (Type I) SN (Aguirre 1999, 
Croft et al. 2000). 

As a final remark, we have seen that SNe with mass 
above M — 15Mq predominantly form amorphous carbon 
grains (the ejecta of these stars have higher C/O ratios); in 
doing so, they use up virtually all the available carbon yield 
in the ejecta. This implies that carbon in the IGM at high 
redshift will be strongly depleted above redshifts at which 
only Type II SNe contribute to the metal enrichment of the 
universe. It will be then posible to test this prediction once 
a sample of target sources at z ^, 5 will become available for 
absorption line studies. 



gas phase mainly via the channel H + e~ — > H~ + hv, at 
rate kg (the rate coefficient kg is given in Abel et al. 1997); 
formation via the Hj channel, when included, is found to be 
negligible in our case. Therefore the formation rate in the 
gas phase is 1Z ~ kgn H -nu. The formation rate on grain 
surfaces is instead given by IZd — 0.5(7C s cr)ndn.sr, where 7 
is the sticking coefficient, c s is the sound speed in the gas, 
and a is the grain cross section. The equality between the 
two rates can be cast into the following form: T> = O.lVTxe, 
where T) is the dust-to-gas ratio normalized to its Galactic 
value, T is the gas temperature and x £ the gas ionization 
fraction. For typical parameters of the PopIII objects (Cia- 
rdi et al. 2000), H2 production on dust grains becomes dom- 
inant once T> is larger than 5% of the local value. With the 
dust yields calculated above, we then conclude that only 
about 50 SN are required to enrich in dust to this level a 
primordial object. Clearly, early dust formation might play 
a role in the formation of the first generation of objects. 

Even in larger galaxies, which will form later on when 
the overall metallicity and dust levels in the universe have 
increased, dust will be at least as important. For example, 
on the scale of the molecular clouds in a galaxy, it will pro- 
vide the opacity to stop the infall on forming protostars, 
hence possibly changing the properties of the IMF, and to 
allow the cloud to self-shield from damaging H2 photodisso- 
ciating UV radiation. Finally, dust photoelectric heating is 
known to be the major heating source for the diffuse ISM, 
and hence partecipating to the onset of its observed multi- 
phase structure (Ricotti, Ferrara & Miniati 1997, Spaans & 
Norman 1997). 



This work was completed as one of us (AF) was a Vis- 
iting Professor at the Center for Computational Physics, 
Tsukuba University, whose support is gratefully acknowl- 
edged. 
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